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Gravitational Microlensing Near Caustics I: Folds 

B. Scott Gaudi 1,2 and A. O. Petters 3,4 

ABSTRACT 

We study the local behavior of gravitational lensing near fold catastrophes. 

Using a generic form for the lensing map near a fold, we determine the observable 
properties of the lensecl images, focusing on the case when the individual images 
are unresolved, i.e., microlensing. Allowing for images not associated with the 
fold, we derive analytic expressions for the photometric and astrometric behavior 
near a generic fold caustic. We show how this form reduces to the more familiar 
linear caustic, which lenses a nearby source into two images which have equal 
magnification, opposite parity, and are equidistant from the critical curve. In 
this case, the simplicity and high degree of symmetry allows for the derivation of 
semi-analytic expressions for the photometric and astrometric deviations in the 
presence of finite sources with arbitrary surface brightness profiles. We use our 
results to derive some basic properties of astrometric microlensing near folds, in 
particular we predict for finite sources with uniform and limb darkening profiles, 
the detailed shape of the astrometric curve as the source crosses a fold. We 
find that the astrometric effects of limb darkening will be difficult to detect with 
the currently planned accuracy of the Space Interferometry Mission (SIM) for 
Galactic bulge sources; however, this also implies that astrometric measurements 
of other parameters, such as the size of the source, should not be compromised 
by an unknown amount of limb darkening. We verify our results by numerically 
calculating the expected astrometric shift for the photometrically well-covered 
Galactic binary lensing event OGLE-1999-BUL-23, finding excellent agreement 
with our analytic expressions. Our results can be applied to any lensing system 
with fold caustics, including Galactic binary lenses and quasar microlensing. 
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Subject headings: astrometry—stars: binaries, fundamental parameters—gravitational 
lensing 


1. Introduction 


Gravitational lensing has proven to be an exceptional tool for studying a diverse set of 
astrophysical phenomena. Its utility is due, at least in part, to the fact it operates in a number 
of qualitatively different regimes. The term strong lensing, or macrolensing, is usually applied 
when a distant source (typically cosmological) is lensed into multiple, resolved images by an 
intervening mass, such as a foreground cluster or a galaxy. Weak lensing is used to refer to 
the case when multiple images are not created, and the gravitational field of the intervening 
matter serves only to slightly distort the image of the source. For most applications of both 
strong and weak lensing, the source, lens, and observer can be considered static. The term 
microlensing is often used to describe the case when multiple images are created, but are not 
resolved. Typically the separation of the images created by a gravitational microlens are of 
order the Einstein ring radius, 
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Here M is the mass of the lens, D is defined by, D = D 0S D 0 i/Di s , and D os , D a i, and D[ s 
are the distances between the observer-source, observer-lens, and lens-source, respectively. 
In cosmological contexts, angular diameter distances should be used. When d E is less than 
the resolution, individual images in general cannot be distinguished. Due to the small scale 
of 9e, it is typically not a good approximation to assume that the source, lens, and observer 
are static. Therefore the lensing properties can be expected to change on timescales of order 
the Einstein ring crossing time, 
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where v± is the transverse speed of the lens relative to the observer-source line-of-sight. The 
standard observables in gravitational microlensing are therefore the time rate of change of 
the total magnification and center-of-light (centroid) of all the microimages. There are two 
different regimes where microlensing has been discussed: quasar microlensing (Wambsganss 
2001) and microlensing in the Local Group (Paczynski 1996). 


In the Local Group, gravitational microlensing occurs whenever a massive, compact 
object passes close to our line of sight to a more distant star. Microlensing was originally 
suggested as a method to detect baryonic dark matter in the halo of our galaxy (Paczynski 
1986), but has been developed and applied as an important tool in studying a number of 
astrophysical topics, including the stellar mass function (Gould 1996), extrasolar planets 



(Mao & Paczynski 1991), stellar atmospheres (Gould 2001), and stellar multiplicity (Alcock 
et al. 2000; Udalski et al. 2000). The only microlensing effect currently observable is the 
magnification of the background source as function of time. This is because, for typical 
distances in the Local Group, the angular Einstein ring radius is 8-^ ~ 1 mas( M/M & ) 1 / 2 , and 
therefore too small to be resolved with current instruments. The timescale for a microlensing 
event is t^ ~ 40 days. In general, it is much easier to determine the center-of-light of an 
image than it is to resolve it. Thus future interferometers, such as the Space Interferometry 
Mission (SIM), although still not able to resolve separations of (P(mas), should be able to 
measure the centroid of all the images to much better than this, perhaps even down to 10 /zas 
in the case of SIM. Such accuracy is sufficient to easily detect the motion of the centroid of 
the images created in a microlensing event, which is also of order 0 e- This regime is typically 
referred to as astrometric microlensing, as opposed to photometric microlensing when only 
the total magnification is observable. 

Astrometric microlensing has a number of important applications. By combining ground- 
based photometry of microlensing events with photometry and astrometry from an astro- 
metric satellite on an Earth-like orbit, the masses of microlenses can routinely be measured 
(Paczynski 1998; Boden, Shao & van Buren 1998; Gould & Salim 1999), allowing for the 
determination of the compact object mass function in the bulge, including stellar remnants 
(Gould 2000). Astrometric information alone allows for the precise (few %) measurement 
of the masses of nearby stars (Paczynski 1995). Finally, for a subset of events, it will be 
possible to obtain precision measurements of angular diameters of stars in the Galactic bulge 
using astrometric information (Paczynski 1998; Gaudi, Graff, & Han 2002). 

Binary microlenses have proven to be enormously useful in photometric microlensing 
studies. This is primarily because binary lenses exhibit caustics: closed curves on which 
the mapping from the lens plane to light source plane becomes critical, and the point- 
source magnification becomes formally infinite. Regions near caustics exhibit large, rapidly 
changing (with respect to source position) magnification, and are therefore useful not only 
for providing a large source flux, but also high angular resolution. However, binary lenses 
have also proven to be difficult to study both theoretically and observationally. This is 
partially because the lens equation, which describes the mapping from the lens plane to the 
light source plane, is equivalent to a fifth-order complex polynomial in the source position 
(Witt 1990), and therefore is not analytically solvable in general. Furthermore, care must be 
taken when considering finite source effects near caustics due to the divergent magnification. 
However, considerable progress can be made when one realizes that the smooth arcs (away 
from cusps) of caustics that arise in nearly equal-mass binary lenses are well approximated as 
simple linear fold catastrophes, which have generic, simple, and most importantly, analytic 
behavior. Thus the caustics of binary lenses can be analyzed analytically or semi-analytically 
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without reference to the global (and non-analytic) topology of the general binary lens. In 
particular, a simple equation for the magnification of a source near a fold exists (Schneider et 
ah 1992), which has been used in a number of important applications including binary-lens 
fitting (Albrow et ah 1999), stellar atmospheres (Gaudi & Gould 1999), and caustic-crossing 
predictions (Jaroszynski & Mao 2001). 

In contrast to the astrometric behavior of single lenses, which is analytic and has been 
quite well studied (Walker 1995; Jeong, Han & Park 1999), there have been only a few 
preliminary studies of the astrometric properties of binary gravitational lenses (Han, Chun, 
& Chang 1999; Chang & Han 1999; Gould & Han 2000). It is known that astrometric 
binary-lens curves exhibit complex behavior, including instantaneous O(0e) jumps in the 
image centroid trajectory that occur when a point source crosses a binary-lens caustic and 
two highly-magnified images appear in a position unrelated to the position of the centroid 
of the other three binary-lens images. The generic behavior of these centroid jumps, or how 
they are altered by finite source effects, is not understood. As is the case for photometric 
microlcnsing, the astrometric behavior of binary lenses will likely prove quite useful for 
several applications. The usefulness of binary lenses is primarily related to the complex image 
centroid trajectories and the large centroid jumps. Although, in general, these properties 
do not allow one to measure any additional parameters over the single lens case; they do 
allow one to measure these parameters much more easily. In particular, Graff & Gould 
(2002) have shown that lens mass measurements can be made to a given accuracy with 
1-3 orders of magnitude fewer photons with caustic-crossing binary-lens events than with 
single lens events, thus greatly reducing the resources required to achieve one of the primary 
proposed science goals of astrometric microlensing. Caustic crossing binary-lens events are 
also enormously useful for measuring the angular radii of microlensing source stars in the 
bulge, for two reasons. First, the expected ratio of binary-to-single lens events for which the 
source star is resolved is a factor of > 4 for giant sources and > 10 for main sequence sources. 
Furthermore, the large and complex centroid shifts expected for caustic-crossing binary-lens 
events makes the requisite astrometric measurements easier. Gaudi, Graff, & Han (2002) 
have shown that, by combining accurate ground-based photometry with a handful of precise 
astrometric measurements, caustic-crossing binary-lens events should yield ~ 5% stellar 
radius measurements with reasonable expenditure of resources. 5 Thus, given the importance 
of caustic crossing binary-lens events, an analytic study of the generic behavior of astrometric 
microlensing near folds would prove quite useful. 

In quasar microlensing, the separate macroimages of a quasar that is multiply-imaged 


5 Although it is possible to measure stellar angular sizes using astrometric information alone, this generally 
requires very densely-sampled measurements, since the source is resolved only for a short time. 
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by a intervening galaxy or cluster also feel the combined, non-linear effect of individual 
point masses (i.e. stars) in the macrolensing object that are near the macroimage position. 
The individual macroimages are in fact composed of many, unresolved microimages with 
separations of order the Einstein radius of a M ~ M & object at cosmological distances, 
6 *e — l/iasfM/Mg) 1 / 2 . The typical timescale for the source to cross an angle of 9e is — 
15 years; however microlensing light curves should show structure on much smaller timescales 
due to the combined effects of many individual microlenses. Since it was first discussed 
by Chang & Refsdal (1979), cosmological microlensing has been studied theoretically by 
numerous authors (see Wambsganss 2001 and references therein), and detected in at least 
two systems (Q2237+0305, Irwin et al. 1989; Corrigan et al. 1991; Wozniak et al. 2000, 
B1600+434, Koopmans & de Bruyn 2000). Observations have been used to place constraints 
on, e.g., the size of the emitting region of quasars (Wambsganss, Schneider, & Paczynski 1990; 
Wyithe, Webster, Turner, & Mortlock 2000) and the mass function of microlenses (Schmidt 
& Wambsganss 1998; Wyithe, Webster, & Turner 2000; Koopmans & de Bruyn 2000). 

Quasar microlensing differs markedly from microlensing in the Local Group in that the 
surface mass density in units of the critical density for lcnsing (the “optical depth 1 ') is of 
order unity, rather than 0( 10 -6 ) for the Local Group. In the high optical depth regime 
the lcnsing effects of the individual microlenses add nonlinearly, resulting in a complex 
caustic network. Due to this nonlinear behavior and the large number of lenses typically 
involved, calculation of the observable properties of such a lensing system is difficult and time 
consuming. Although in the high optical depth regime the caustics often exhibit considerably 
more complicated global behavior than the caustics of binary lenses in the Local Group, 
it is still the case that the smooth arcs (away from cusps) of the caustics are locally well- 
approximated by generic fold catastrophes. This fact, combined with a simple formula for the 
magnification near folds, has been exploited by numerous authors to quickly and efficiently 
calculate various observable properties of quasar microlensing (Wambsganss & Paczynski 
1991; Lewis & Ibata 1998; Wyithe & Webster 1999; Fluke & Webster 1999). 

The observable effects of quasar microlensing have been limited to the relative magnifi¬ 
cations of the various macroimages as a function of time. As with Local Group microlensing, 
astrometric effects should also be present. The centroid of the individual macroimages should 
vary as a function of time, particularly when new images are created or destroyed when the 
source crosses a caustic. This effect has been studied by Williams & Saha (1995) and Lewis 
& Ibata (1998). In particular, Lewis & Ibata (1998) predict that magnitude of the centroid 
shift for the Q2237+0305 system (apparent magnitude R < 18.5) can be as large as ~ 5Cffias, 
and thus potentially observable with SIM. They also note that the magnitude of the cen¬ 
troid shift is often correlated with the magnitude of the change in total magnification. The 
analytic results presented in this paper may prove useful for this application. 
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Here we study the generic, local behavior of microlensing near fold catastrophes. In 
§2 we present an analytic study of the photometric and astrometric behavior near folds. 
We begin with the equations that describe the mapping near a fold caustic in §2.1, and 
use these to derive the behavior of a point source near a fold. We extend this analysis to 
finite sources in §2.2, and limb-darkened sources in §2.3. In §2.4 we show how and when our 
generic parabolic fold form reduces to the more familiar linear caustic. In §2.5 we use our 
analytic results to derive some generic results about the astrometric behavior near folds. We 
verify the applicability of our results in §3 by numerically calculating the photometric and 
astrometric behavior of one well-observed binary-lens event. We find excellent agreement 
with our analytic formulae. Finally, we summarize and conclude in §4. 

Our goal is to provide a thorough, comprehensive study of gravitational microlensing 
near fold caustics. Although our study is interesting in its own right, the primary utility 
of the results presented here is their potential application to the topics mentioned in the 
previous paragraphs. A prescription for how specifically our results can be applied to these 
topics is beyond the scope of this paper, but we will make general comments along these lines 
over the course of the paper, and more specific comments in §2.5. We are currently preparing 
a complementary, similarly detailed study of microlensing near cusps. Combined with this 
study, we will have a reasonably thorough and complete understanding of the local behavior 
of microlensing observables near all stable gravitational lcnsing singularities. We note that 
some of the results derived here, particularly the results on the photometric behavior near 
folds, have been presented elsewhere (see, e.g. Schneider et ah 1992, Petters et al. 2001, and 
Fluke & Webster 1999). We include those results here for the sake of completeness. 


2. Analytic Considerations 
2.1. Lensing Near Fold Caustics 

2.1.1. Global Lensing Map 

For a general gravitational lens, the lensed images 9 of a source u, are given by the 
solutions of the lens equation, which is the mapping u —> 6: 

u = rj(9) = 9 — cx(9), (3) 

where a = V"0, and tfj is the projected Newtonian potential of the lens, 

iP( 0 ) = - [ d9'n{9') In \9 - 9' 

11 Jr? 


( 4 ) 
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Where k = E(0)/E cr , E(0) is the surface density of the lens, and 


. g2 D os 
47 tG DoiDis 


( 5 ) 


is the critical surface density for lensing. Note that V 2 ^ = 2k. We are assuming that 
0 = r j ( 9eD 0 i ) and u = s/(9 eD os ), where r and s are the proper vector positions in the lens 
and light source planes, respectively. The mapping r] can produce multiple images of the 
source u. The magnification of lensed image 0, is 

" (9i) = (6) 


where Ajj is the Jacobian matrix of the lensing map r/. The critical curve is the set of 
positions 0 C such that the determinant of the Jacobian matrix vanishes, i.e. where 


J = det[Ar,(6 c )]=0, (7) 

and the caustics are u c = rj(6 c ). In the case of microlensing, the individual images are by 
definition unresolved, and thus it is useful to define the total magnification 

Ikot ^ ^ A i 
i 

where the sum is over all images. The center-of-light, or centroid O c i{u) of the images is 
simply the magnification weighted sum of the image positions, 

l^tot 

For simplicity, we will focus primarily on the quantity 6d. However, it is important to note 
that, in general, centroid measurements will be made with respect to the unlensed source 
position, thus the observable is 


Od(u) 


(9) 


( 8 ) 


60 d (u) = 0 d (u) - u. ( 10 ) 

Also, while the angular variables we will be working with will be in units of 0e, the astrometric 
observables are in physical units, such as arcseconds. To convert to observable quantities, 
all angular quantities must be multiplied by 9e■ For example, the observable centroid shift 
is given by 5ip cl = 9 E S0 d . 
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2.1.2. Lensing Map Near Folds 

We now derive the generic behavior of the photometric and astrometric properties of 
gravitational lensing near folds. We will present onr derivations in some detail, in order 
to document the approximations and simplifying assumptions that are implicit in the final 
analytic expressions. In Figure 1, we provide an illustrative example of the basic properties of 
lensing near a fold. We will refer to this figure repeatedly during the course of the derivations. 

Suppose that the lensing map rj sends the origin to itself (which can always be accom¬ 
plished by appropriate translations) and a fold caustic curve passes through the origin. By 
Taylor expanding the gravitational potential 0 about the origin, one can find an orthogonal 
change of coordinates that is the same in the lens and light source planes such that the 
lensing map 77 can be approximated by the following mapping in a neighborhood of the 
origin (Petters, Levine, & Wambsganss 2001, pp. 341-353; Schneider, Ehlcrs, & Falco 1992, 
p. 187): 

= ad 1 + — d\ + cOid- 2 , U 2 = —Of + bd 162 + —$2 > ( 11 ) 

where (di,d 2 ) and (ui,u 2 ) denote the respective coordinates in the lens and light source 
planes, and 


a = 1 - 0n(O) 7 ^ 0, 6 =-■ 0122 ( 0 ), c=-0n 2 (O), d = - 0222 ( 0 ) 0 0. (12) 


Here the subscripts refer to the partial derivatives of 0 with respect to the original global 
Cartesian coordinates of the lensing map. For the example in Figure 1, we have adopted 
a = 5, b = 1, and c = — d = —0.5. 


The Jacobian matrix of equation (11) is 


A = 


a T cd 2 cd\ T bd 2 
cd\ + bd 2 bd] + dd 2 


(13) 


The critical curve is given by 


J = det A = (a + cd 2 ){bd\ + d# 2 ) — (cd± + bd 2 ) 2 = 0. (14) 

The tangent line to the critical curve at the origin is given by 

0 = 6 ■ VJ(0) = abd\ + add 2 , (15) 


that is, 


02 =—A 
d 


( 16 ) 
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since a ^ 0 and d ^ 0. Substituting 9 2 = —b6i/d into equation (11) yields 

ui = a6>i + ^ Q& 2 - cd'j 9\ ~ a6>i, w 2 =-^(cd - 6 2 ) d 2 . (17) 

Note that in the expression for iq, the term 6 1 ! dominates 6^ near the origin. Inserting 
9\ = Ui/a into u 2 above, we see that the tangent line at the origin of the critical curve is 
mapped into a parabola (Schneider et al. 1992; Fluke & Webster 1999): 

C(u\, u 2 ) = 2a 2 d{u 2 — eu\) = 0. (18) 


where we have introduced the combination of local derivatives of i/j: 

_ cd — b 2 
2a 2 d 


(19) 


We show in §2.4 that \e\ is one-half the curvature of the caustic at the origin. Thus when 
\e\ -C 1, the caustic can be approximated as 2 a 2 du 2 = 0, i.e. the tq-axis. For the example 
shown in Figure 1, |e| =0.05. 

Since the tangent line (Eq. 16) approximates the critical curve near the origin, the 
parabola (Eq. 18) approximates the caustic near the origin. See Figure l(a,b). Now multi¬ 
plying u 2 in equation (11) by 2 a 2 and substituting 9i = u\/a yields 


C(ui, u 2 ) — (b u\ + ad 0 2 ) 2 . 


( 20 ) 


Hence, if a light source is located at a position (wi, - ^) where C(ui,u 2 ) < 0, then there is no 
lensed image locally, while a source with C(ui,u 2 ) > 0 has at least one image. 

For an n-point mass lens, we have c = —d. Consequently, if d > 0, the parabolic caustic 
lies in the lower-half plane locally. In addition, the region above the parabola is such that 
sources lying there have double images locally (Fig. la). If d < 0, then parabola is in the 
upper-half plane with the region below the parabola yielding double images locally. In other 
words, the caustic is locally convex (see Petters et al. 2001, Sec 9.3 for a detailed treatment). 


2.1.3. Image Positions of Sources Near Folds 


Let us determine the images for C(ui, u 2 ) > 0. Equation (20) is equivalent to a quadratic 
equation in 0 2 , 


0 = [a 2 d 2 ] 6 2 + [2abdui] 9 2 + \b 2 u\ — C(ui,u 2 )}. 


( 21 ) 
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The solutions are 6 


e 2 = 


-bUi ± y/C(ui,U 2 ) 


ad 


The expression for U\ in equation (11) yields 

1 


6 1 = 


a T cd 2 

Ignoring terms of order 3 or higher, we obtain 


«i - - 2 <>l 


( 22 ) 


(23) 


«i = - - 

a cr 2a 


(24) 


where the approximation (a + c# 2 ) 1 — a *(1 — 062 / 0 ) was employed. Substituting equa¬ 
tion (22) into (24) and keeping only terms that are linear in u\ and u 2 , it follows that 


0 i = 


d/a 1 — bu 2 
ad 


Hence, a source with C(ui,u 2 ) > 0 has two (opposite parity) images given by 

o± = (^±,1.^,2) - 2 ^ ^ dui - bu 2 , -bui ± v/eTyytjj. 


(25) 


(26) 


Figure l(a,b) illustrates the mapping from source to images near a fold. 

Consider a point w c on the caustic, i.e., C(u c ) = 0. Locally there is one image of u c , 
located on the critical curve 7 at the position 

1 

Oc = —;(du 1 - bu C)2 , -bu cl ). (27) 

ad 


Now consider a source at it = (ui,u 2 ) inside the caustic (i.e., C(u) > 0) and let u* = 
(uh,u* 2 ) be the point on the caustic with the same horizontal position as u , i.e., u* c l = U\ 
and u* c 2 ^ u 2 . The vertical separation of the source from the caustic is then the difference 
between the '(^-coordinates of u and u*: 


Aw 2 ,_l = u 2 - e(u* c l ) 2 


C(ui,u 2 ) 
2 a 2 d 


(28) 


6 Note that we do not need to include sign (ad) in front of the square root since ±sign(ad) = ±1. 

'That the image position in equation (27) corresponds to a point on the critical curve can be seen by 
noting that on the caustic w c> 2 = eu 2 c l5 and close to the origin, U\ dominates u\. and thus the image position 
is 0 C = (dw c 1 , — bu c ,i)/{ad), which is indeed a point on the critical line (see Eq. 16). 



11 


Note that Am 2i i > 0 if and only if d > 0 (since the source is located where C(ui,u 2 ) > 0). 
Using equations (26) and (27), we find that the distances of the images 0±(ui,u 2 ) from the 
point 6 C on the critical line are the same, 


A6» 2 ,_l = |0±(tti,tt 2 ) - 0 C | = 


2A u 


2,_L 


d 


1 + 


b 2 

2 a 2 d 


1 1/2 


A u 


2 ,X 


(29) 


Near the origin, the first term in brackets dominates over the second, and thus, 

A$2,± = 


2A u 


2,-L 


d 


(30) 


2.1.4- Magnification of Sources Near Folds 


Using the earlier tangent line approximation to the critical curve at the origin, we Taylor 
expand the Jacobian determinant J about the origin to first order: 


J{9i,9 2 ) — Ji9i(O)0i + </i9 2 (O)02 — ab 9\ + ad 9 2 , (31) 


which is equal to the right-hand-side of (15). Inserting the expressions for 9\ and 9 2 from 
equation (26), we obtain 

J(9i, 9 2 ) = — —w 2 ± y/ciu^u^}. (32) 

Keeping only the lowest order terms in u j and u 2 gives, 

\J{0±)\ = VC^ufi. (33) 


In terms of the vertical separation between the source and the caustic (Eq. 28), the magni¬ 
fication of the images are given locally by 


= M#±) 


l 

wai 


1 n*7 

2 V A u 2y 


(34) 


where Uf — 2 /(a 2 d). Thus the two images have the same magnification if the source is 
sufficiently close to the caustic. The total magnification Uf = 2/i ± of the two images is 
simply 


jif 


/ 


u f 

A u 2 ,± ’ 


(35) 


which agrees with the expressions derived by both Schneider, Ehlers & Falco (1992, pl90) 
and Fluke & Webster (1999). Equation (35) is the well-known result that the magnification 
varies inversely as the square-root of the distance of the source from the caustic. It is not 
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often appreciated however, that this property holds for the general parabolic fold form as 
well, provided that the vertical distance is used, rather than, for example, the minimum 
distance between the source and the fold caustic. Rigorously, the vertical distance is the 
distance between the source and the caustic in the direction perpendicular to the tangent 
line of the caustic at the origin, where the origin is defined as the point around which the 
potential ip is Taylor expanded. In the limit of a straight fold (|e| <C 1), the vertical distance 
and the minimum distance are equivalent. 


2.1.5. Image Centroid of Sources Near Folds 


Write the source position u in terms of the fold local coordinates, i.e., u = (u 1 ,u 2 )- 
There are no local images for source a position with C(u) < 0. From equations (9) and (26), 
and since /a + = for sources inside the caustic the centroid is given locally as 


o.f( u ) = 7} ( 0 + + e ~) = ~ ~ bu ^ 


for C{u) > 0. 


(36) 


We express the rectilinear motion of the source as follows: 

u{t) — u c + (t — t c )u, (37) 

where u c is the position at which the source intersects the caustic at time t = t c and ii is 
the constant angular velocity vector of the source, 

" cos (p sin <p ' 


(38) 

Recall that £e = $e D a i/v±, where v_i is the transverse speed of the lens relative to the 
observer-source line-of-sight. Here (p is the angle of the source’s trajectory with respect to 
the Mi-axis. Note that the iq-axis does not necessarily coincide with the caustic, and thus 
<p> is not necessarily the angle of the trajectory with respect to the tangent to the caustic at 

U r - 

We shall assume that t > t c corresponds to the source’s trajectory lying in the double- 
image region. In that region, the centroid follows a straight line locally: 

df(t) = 0/, c + (t - £ c )0(/, c ), for t>t c (39) 

where 0f tC = 0f(u c ) and 0 (/, c ) = d0f/dt\ t=tc = 0f{u) (both constant vectors). Note that 
since the slope of the centroid line Of is, tan (p c = (—b cos <p)/(d cos <p — b sin d), where <p c is 
the angle between the centroid line and the di-axis, it follows that 0 (f, c ) can be expressed as 

-bcoscp 




-(cot (pc , 1). 


ad £e 


(40) 
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2.1.6. Global Magnification and Centroid 


Fold caustics do not, of course, exist in isolation. They are tied to the global properties 
of the lens in consideration. For practical purposes, we therefore consider images created by 
the lens that are not associated with the fold under consideration. We define /x 0 to be the 
total magnification of all the images not associated with the fold, and we define 6 0 as the 
centroid of all these images. We will assume that there is locally only one fold caustic, and 
that all the other image magnifications and positions are only slowly varying functions of 
the source position. The total magnification is then 


dtot, df + do 



@(Aw 2 ,_l) 


+ do- 


(41) 


See Figure 1(c). Here 0 is the Heaviside unit step function (i.e., 0(x) = 1 for x > 0 and 
0(x) = 0 if x < 0); it accounts for the fact that sources below the caustic (with A-u 2 ,_l < 0) 
have no images locally (i.e., near the critical curve). Since fi + = d- = \Xf/2 , we have, 

@ci — - \df@f + do@o] (42) 

dtot 


where 


0/H 


l 

ad 


(du\ — bu- 2 , —6mi)0(Am 2i j_). 


(43) 


We can calculate the dependence of the observables dtot and 0 c i on time by assuming 
a rectilinear source trajectory and replacing u in equations (41) and (42) by equation (37). 
We first Taylor expand 9 0 and /i 0 about the time of the caustic crossing, t = t c , keeping 
terms of first order in t — t c : 


Qq — $o,c T (t — £ c )0(o,c)) 
do do ,c T (t t c )/i(o,cp 


(44) 

(45) 


where 9(o,c) = d0 0 /dt\ t=tc and fi(o,c) = d/x 0 /df| t=tc . From equations (28) and (37), Am 2 j_ is 
given by 

t — t r 


A« 2 ,i = 


tv, 


■ sin <f) 


1 — e cot <f> I 2u C; i H-- cos <f> 


t-R 


(46) 


The magnification as a function of time is then, 

1/2 


dtotif ) 


t 


f 


t — t f 


1 — e cot <t> ( 2 u c i H- -u f cot <f 


- 1/2 


<d(t-t c ) 

+ dO,c + d(0,c)(t — t c )i 


( 47 ) 
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where tf = (uft E )/ sin (j) is the effective rise time of the caustic crossing, and Uf — 2 /(a 2 d) 
defines the characteristic rise length as before. Notice that when the curvature of the caustic 
is small, i.e., when |e| <C 1, the magnification associated with the fold reduces to the more 
familiar form /.if = [{t — t c )/tj]^ 1 / 2 @(t — t c ). The centroid of all images is 

@cl(t) — ——' | — t c ) + 00^0,c] + [/hf0(/,c)©(^ — tc) + ^0^0 ,c](t — t c )| , (48) 

where 6f c and #(/, c ) are defined below equation (39). Figure l(e,f) illustrates the behavior 
of the two components of 6 c i{t) as a function of time, whereas Figure 1(g) shows 9d ,i versus 

9 cl,2 ■ 


2.2. Finite Sources 


2.2.1. Finite Source Magnification 


The results of the previous section assumed a pointlike source. This results in an 
astrometric curve that exhibits an instantaneous jump from 6 o >c to 0f tC at t — t c . All real 
sources will have a finite extent which will smooth out the discontinuous jump. For a finite 
source, the magnification is the surface brightness weighted magnification integrated over 
the area of the source, 


f s = fp d uS(u)/i(u) 
f D duS(u) 


(49) 


where D is the disc-shaped region of the source and S(u) is the surface brightness of the 
source. Let S be the average surface brightness of the source, S = (npl)^ 1 f D d uS(u). Here 
p * = 0*/9 e , is the angular source radius 9 * in units of 9 E . The denominator in equation (49) 
is then simply tt p*S. Dehne Sn = S(u)/S to be the normalized surface brightness. Also 
define a new set of source plane coordinates such that 


u 


/ 


u — u 


cn 




(50) 


where u cn is the position of the center of the source. Then for a point u inside the disc 
source H, we have \u'\ < 1. Equation (49) becomes 



(51) 


All of the preceding results apply to generic parabolic fold catastrophes. However, in 
order to continue making significant progress analytically, we must make the following simpli¬ 
fying assumption: We furthermore assume that \e\ -C 1, and thus the caustic coincides with 
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the u i-axis and Au 2 ^± = U' 2 ■ This considerably simplifies the form for the fold magnification 
Pf- 


Let z be such that zp * is the vertical separation of the center u cn of the source from 
the caustic, i.e., u CH:2 = zp *. The source D is in the upper-half plane (i.e., u 2 > 0) if and 
only if z > 1. If z — 1, then D just touches the caustic, while for 0 < z < 1 a portion of 
D is below the caustic with the center of D on the caustic for z = 0. The center of D lies 
below the caustic for z < 0, with a portion of D still above the caustic for —1 < z < 0 and 
D completely below the caustic for z < — 1. Since the fold magnification pf is nonzero only 
for points u in D that he in the upper-half plane, we get 

, fu f \ 1/2 0(1 + 2) . , 


04 + z 

By equations (51) and (52), we hnd for a fold caustic and arbitrary surface brightness profile, 

0(1 -f z) fV 1 -^ 

ClUo , = 


fs j 

0 1 


z = 


= (?l) 
P*J 


1/2 


1 

7T 


V U 2 + Z 00+^1 


.d u[ Sjv(ui,i4) 


(53) 


• max(- 2 ,-l) 

For a uniform source, i.e., Sn(u') = 1, this simplifies to (Schneider et al. 1992; Albrow et al. 
1999), 

= fe) 1/2 

,P*J 


where 


G r 


(Z) = 7T- 1 / 2 - 


p f [Z) = 

(n + 1)! 


Go(z), 


(54) 


(1 


, dx0 — X - . 0(1 + z). (55) 

(n + 1/2)! y max( _ z ,_i) (x + z) V 2 

Note that G 0 can be expressed as an elliptic integral. Figure 2 shows G 0 (z) for —2 < 2 < 2. 

For small source sizes 0* <C &e, the magnification p 0 of the images not associated with fold is 

a slowly varying function of u over the source, and thus 7 r” 1 f D clu'Sjv(u')po(u') = po(u cn ). 

Therefore, the total finite source magnification is just 


2\n+l/2 


Ptot ~ P U f + P'0,cm (56) 

where p 0tCn = p 0 (u cn ). In analogy to the point source case (§2.1), we can expand p 0 about 
the time of the second caustic crossing. For a source with rectilinear motion, the total finite 
source magnification in the time domain is then, 

Ptot(t) = (a/) ^0 ( + Po,c +P(o,c)(t ~ t c ), (57) 

where tf = (0)/ sin 0 (effective rise time of the caustic crossing) and At = (p*t e)/ sin (f) is 
the time scale of the caustic crossing, i.e., the time between when the source first touches the 
caustic and when it straddles the caustic. Note that tf/At = Uf/p *. Figure 1(d) illustrates 
the behavior of p// t (t). 
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2.2.2. Finite Source Image Centroid 


For an extended source, the image centroid is the position of the images weighted by 
both the surface brightness and magnification, integrated over the area of the source: 

nfs _ Ei Id d uS(u)0i(u)fj.i(u) 

Zif D d nS(u)Mu) ’ 1 J 

where the sum is over all the microimages. The denominator is simply the total flux, STTpl/ij^. 
Again, for -C $e, the centroid 9 0 (u) of the images not associated with the fold varies slowly 
over the source. This yields ti 1 f D du / Sjv(u / )p 0 ( u ')9o(u') = So( u cn)9 0 (u cn ) (since po(u) is 
also slowly varying over the source). Equation (58) then separates into two terms: 


nfs 

“d 


fs 

fs 
Stot 
fs 

fs 

Slot 


TTplSfL f / 
fs , So,ci 


d uS(u) 9f{u ) pf( 


u 


’ D 


+ 


So ,c 


fs 

Slot 


0 


0 ,cn 




fs 

Slot 


<9 


0 ,cn? 


(59) 

(60) 


where 6 0jCn = 9 0 (u cn ). The hrst term in (60) is the contribution from the two images 
associated with the fold caustic, while the second term is the contribution from all unrelated 
images. For convenience, we defined Of to be the factor within the brackets in equation (59). 


The term 9^ s will now be evaluated. Since u = p*u' + u cn , we obtain 


9f{ y 'd') 9f(u cn ) ~h p*9f{u ). 


(61) 


Define 9f :Cn = 9f(u cn ) 0(1 + z), so 9f :Cn vanishes when the disc source lies completely below 
the fold caustic. 8 Using (52), we find a simple formula for the finite source image centroid, 


ej‘(t) = e Lcn + FFp- e,(H s „,ic Sn ), 

d 

where we have defined 


HsM 

1 

[ du' 2 

0(1 + Z) 1 

du\ u\ 

Sn{ u 1, u 2 

7T , 

J max(— 2 :, —1) 

VU 2 + Z J- 

-\J !-( u 2 ) 2 


1 

[ du’ 2 

u > 0(1 + A> 

/V 1_ 64) 2 

/ du) 

S N (u[,u 2 

7T . 

J max(- 2 ,-l) 

sju 2 + z 



(62) 


(63) 

(64) 


technically, the fold produces no images locally of u cn when z < 0, and thus 0f(u cn ) is not defined for 
z < 0. This apparent discrepancy can be alleviated by simply assuming that 0 /(ui, U 2 ) is a function that is 
defined for all iti,W2, he., 0 f{u\,U 2 ) = ( du\ — bu 2 ,—bu\)/(ad) for all ui,U 2 - 
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In the case of a uniform source, we obtain 9 Hs N (z ) = 0 and ICs N = /Co, where 


JCn(z) = 7T 1/2 


(n + 1)! f 1 


("1 _ r 2\n+l/2 

dxx— -—— @(1 + z). 


{n + 1/2)! 7 max (_ 2 _i) (x + z) 1 / 2 


For a uniform source, the finite source centroid of the two images associated with the fold is 
then, 

Of = O f ,cn ~ b^^ lCoiz) i, (66) 

ad Hf 

where i = (1,0). Inserting the definition of /j,y, this can also be written in the alternate 
form, 

nus a _ bp* tCo(.) . , 

“f —Uf,cn _ j n / n 1 - ( 67 ) 


aus a bp*JC 0 (z) 

Vf ~ &f,cn -J r* t x 1 - 

ad Gq{z) 


Figure 2 shows the functions K.q(z) and Go(z), while Figure 3 depicts )Cq(z)/Gq(z). Figure 
1(g) illustrates the behavior of the image centroid (Eq. 60) for a finite source. 


2.3. Limb Darkening 

In this section we consider the effect of non-uniform sources on the magnification and 
centroid shift near folds. The effect of generic surface brightnesses can be evaluated using 
the general integral forms for the magnification (Eq. 53), and the two components of the 
centroid shift (Eqs. 63 and 64). We will concentrate on a specific form for S(u) applicable 
to stellar sources, namely 

Sn(u) = jl - r 1 - \ (l - ^ f J , (68) 

where u cn is the center of the source. Here V is the limb darkening parameter, which may be 
wavelength dependent. This form was originally introduced by Albrow et al. (1999), and it 
has the desirable property that there is no net flux associated with the limb darkening term. 

Inserting this form for the surface brightness into equation (53), we recover the result 
of Albrow et al. (1999) for the limb darkened magnification 

/ \ 1/2 

p l f = PY + r \^~J [^1/2(2) — Gq(z)\, ( 69 ) 

9 In fact, Hs n (z) = 0 for any profile with the symmetry Sn(u[, u' 2 ) = <Sjv(— u' L , u' 2 ). 



18 


where G n (z ) is defined in equation (55), and Gq(z), Gi/ 2 (z), and Gi/ 2 (z) — Go(z) are shown 
in Figure 2. 

Since this form of the surface brightness profile is symmetric, Sn(u[, u' 2 ) = Sn{—u\ , i// 2 ), 
the integral Hs N vanishes. Inserting equation (68) into the general form for fCg N , we find 
ICs N = /Co + r(/Ci/2 - /Co), and thus, 


e 


Id 

f 


e 


f,cn 


b ^/u]~pl 
ad (i l f 


\K 0 (z) + r(/c 1/2 - /C 0 )] 


which can also be written, 


0 


Id 

f 


e 


f, 



~ /C 0 (;g) + r(/C 1/2 -/Cp) ' 

_G 0 (z) + r(Gi/2 — Go). 


(70) 


(71) 


Figure 2 shows the functions /Co(^), ICi/ 2 (z), and ICi/ 2 (z) —ICo(z), while Figure 3 depicts the 
term in brackets above for several values of T. 


2.4. Simple Linear Folds 

The majority of the results presented in §2 (with the exception of the assumption made 
in §§2.2 and 2.3 that the caustic is coincident with the iq-axis) are applicable for the general 
case of a parabolic fold catastrophe. In some cases, however, it is possible to simplify the 
expressions considerably and recover the more familiar linear fold form. 

To see how this limit may be reached, consider the form for the general fold mapping in 
equation (11). In this form, the source and image plane coordinates have all been normalized 
by 6e, and thus are dimensionless, as are the coefficients a,b,c,d (since ijj is dimensionless). 
It is clear that the caustic will (in general) be appreciably curved when one considers order 
unity variations in ui,u 2 (i.e., on absolute angular scales of order 9e)- In other words, in 
the general case one can expect that all of the coefficients (a, b, c, d) to be of the same order 
of magnitude. Now consider variations on some smaller scale 9 e -C 9e- Renormalizing the 
angular source and image plane variables such that u = ( 9 E /9 6 )u and 6 = ( 9E/9 e )G , the fold 
mapping can be recast in the form 

hi = AO i + —^2 + G9\9 2 , u 2 = — 9\ + B9\9 2 + —#2> (72) 

where the relation between the coefficients are A = a, B = b(9 e /9E), C = c^/^e), and 
D = d(9 t /9E)- Therefore, on scales of order 9 e , we can expect A to be larger than B, C 
and D by a factor 9e! 9 <l . In the Galactic microlensing case, for example, one is typically 
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concerned with variations on the scale of the source, i.e., scales of order 9 e rs #*. In this case 

Oe/ 0 e > 100 . 


Let us consider the curvature of the fold caustic on scales of 0(9 e ). In general, for a 
twice continuously differentiable function / = f(x) on an open interval of R, the magnitude 
of the curvature at a point (x 0 , f(x o)) on the graph of / is given by 


h x o)i 


f"(x o) 

[1 + (/'(io)) 2 ] 3/2 ' 


(73) 


In the coordinates of (72), i.e. on scales of 0(9 e ), equation (18) yields that the fold caustic 
can be expressed as a graph (hi, / (hi)), where 



. CD — B 2 
f(u 1 )=eu 1 , e- 2A2D . 

(74) 

By (73), 

- 2 lei 

(75) 


“ [1 + 4e 2 h 2 ] 3 / 2- 

But 

e = hr e ’ h hr 

V / \ / 

(76) 


Hence, the magnitude of the curvature of the fold caustic can be expressed as 


ih«i)i 


[OA 

l 21 

e 


V^E/ 

' [1 + 4e 2 

m 2]3/ 2 


(77) 


Equation (77) yields that, at the origin and on scales comparable to 9 e , the magnitude 
of curvature scales as 9J9e, and is given by \k(0)\ = 2(# e /#E)e. Note that the radius 
of curvature at the origin is |/c(0)| _1 . In Galactic microlcnsing 9 e ~ 9 * and #*/#e <C 1. 
Therefore it is generally the case that \k\ <C 1, i.e., the curvature of the fold caustic is 
negligible on scales of order and the fold can be treated as linear. 


We therefore ask what happens to the expressions derived in the previous sections in 
the limit that a 3> ( b,c,d ). We find that C(ui,u 2 ) — 2 a 2 du 2 , and thus the caustic, which 
is defined by C(ui,u 2 ) = 0, collapses to the rq-axis (linear fold). The critical curve is still 
given by 9 2 = —(b/d)9\. There are images whenever u 2 > 0, which are located at 


0 ± = 



(78) 


/^zb 


1 

2 V u 2 


with magnifications, 


(79) 
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The centroid of the two images therefore takes on the simple form, 


u l 


»,= -,o • 


(80) 


Including the additional images, and assuming a rectilinear trajectory, the total magni¬ 
fication can be written, 


( t f \ 1/2 

Ahot(^) ( - ~ J Q(t t c ) T /io,c T /i-(o,c)(^ t c ), 


The centroid of all the images has the same form as before (Eq. 48), with 

"cos (f) 


0f ’ c ~ ( a ’ 1 ’ 0 


e (f,c) - 


ait 


-,0 . 


(81) 


(82) 


For finite sources, the centroid in (66) reduces to 0j s = 9f^ cn . The finite source centroid 
then becomes 


fs 

nfs _ q fJ- 0,cn Q 

“cl f s V f,cn + t s ,cni 

riot riot 


(83) 


This is the same result as for point sources (see Eq. 42) if the point source magnification is 
replaced with the finite source magnification. 

Equations (78)-(83) should hold whenever (1) a 3> ( b,c,d ), (2) one is not too close to a 
higher-order catastrophe, i.e., a cusp, and (3) there are no other nearby folds. 


2.5. Some Applications 


We can use the results from the previous sections to derive some generic results about 
the astrometric behavior near folds. We first consider the magnitude of the astrometric 
jump when the source crosses a caustic. It is clear that for a point source, the maximum 
centroid shift is A 9j ump = 1 0f, c — 0o.c|, where 9f c is the point on the critical curve where 
the images merge. Although A 9j ump ultimately depends on the source trajectory and the 
topology of the lens, we can typically expect that A 6j ump is O(0e)- For the finite source 
case, the difference between the centroid position just before (z ~ —1) and just after (z ~ 1) 
the caustic crossing is A 9^ s amp ~ (/ij /) Ad lump . For a uniform source, this is 


Ad 


US 

jump 


i + 



Ad 


jump • 


(84) 
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Adopting typical parameters, po,cn ~ 1, Uf ~ 1, and Go ~ 1, we find the fractional change 
from the point source case to be (A— A9j ump ) / A9 jump ~ —3%(p*/lCT 3 ) 1 / 2 . 

We now consider the magnitude of limb darkening effect on the magnification and cen¬ 
troid shift relative to the uniform source case. From equations (56) and (69), the fractional 
difference between the limb-darkened and uniform-source magnification is, 


S/Md = dkg* = r Gl/2 “ Go 


/C 


G(] 



-i 


(85) 


Early in the caustic crossing (z < 1), the factor a/ p*/uf(fi 0)Cn /G 0 ) is small compared to 
unity, and thus 5/xid ~ r(Gi /2 — Go)/Go. The magnitude of (G 1/2 — Go)/Go is < 20% for the 
majority ( z < 0.8) of the caustic crossing. Near the end of the caustic crossing (z ~ 1) the 
term in brackets begins to dominate as Go —> 0, and thus 5/iui r U f/ P*{Gl/2 — G 0 ), 
which goes to zero as z —> 1. Thus for typical values of Uf,p*, and /i 0iCn , the fractional 
difference from a uniform source is < 0.2T for the majority of the caustic crossing. See 
Figure 3. 

The difference in the centroid due to limb darkening is 

= 0“ - 0'S = 5iM d ^ (0, - 0„). ( 86 ) 

Gtot 

Figure 3(c) shows the prefactor h/Pcz/io,en//4of f° r ho,cn = 4 , Uf = 1, and several values of p*. 
To assess the detectability of the deviation of the centroid shift from the uniform source case 
due to limb darkening in Galactic bulge microlensing events, we now make a crude estimate 
for the maximum magnitude of AO 1 // by adopting typical parameters. At t — t c , we have 
that | Of — 0 O | = A 6jump 1 and is maximized. We will assume that A6 Jurnp ~ 9e ~ 300pas. 
We have just argued that S/iid < 0.2r. Near the end of the caustic crossing (where is 
maximized), n\ d ot ~ /i 0 • Inserting these values into equation (86), we find 


AO l /l\max ~ (60r)pas. (Bulge Lenses). 


(87) 


For typical values of T ~ 0.5 in the optical, \A0 l ( /\ max < 30pas, which is only a factor of 
~ 3—6 larger than the sensitivity expected from SIM. Note that this is the maximum centroid 
shift; inspection of Figure 3(c) reveals that \A0 l //\ is considerably smaller than this maximum 
for the majority of the caustic crossing. Therefore, detection of the astrometric effects of 
limb-darkening will likely be challenging, at least for Galactic bulge lensing events, with 
the currently planned accuracy for SIM. On the other hand, the fact that the astrometric 
effects of limb darkening are small implies that limb-darkening can generally be ignored. 
Thus predictions for the signatures of other effects, such as the finite source effect itself, are 
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robust. In other words, measurements generally should not be compromised by an unknown 
amount of limb-darkening. 

Consider two observers that are not spatially coincident. The source trajectories as 
seen by the two observers will be displaced relative to each other by an amount Su, the 
magnitude of which is |<Sit| = a© s /rE, where a© s is the component of the separation between 
the two observers perpendicular to the line-of-sight, and r e = O^D is the Einstein ring radius 
projected to the observer plane. If Su is sufficiently large, then both the photometric and 
astrometric behavior of the event will be measurably different between the two observers, 
an effect commonly known as parallax. Since a® s is known, this observed difference can 
be used to infer Pe, providing additional physical constraints on the properties of the lens. 
For example, when combined with a measurement of 9 e from the astrometric centroid shift 
itself, the mass of the lens can be determined, M = (c 2 /4G)rE#E- For a given |5w|, parallax 
effects are typically largest near caustics, because the magnification and the centroid vary 
rapidly with respect to source position. Thus caustic crossings are ideal for use in measuring 
parallax effects (Hardy & Walker 1995; Gould & Andronov 1999). If |5it| is sufficiently small 
that the behavior near the fold caustic crossing for both observers can be described by the 
same local expansion of the lens mapping, than the results derived here can be used to fit 
the astrometric and photometric behavior for both observers and derive the parallax effects 
without regard to the global behavior of the lens. We note that, for generic fold caustics, 
any significant displacement Su, regardless of the its orientation relative to the caustic, will 
result in a difference in the observable behavior. However, for linear fold caustics (§2.4), 
only displacements perpendicular to the caustic will result in significant differences in the 
magnification and centroid. Thus, for generic linear fold caustics, only a projection of Pe is 
measurable from the local caustic behavior (Graff & Gould 2002). 

During a fold caustic crossing the source is resolved, altering both the photometric and 
astrometric behavior with respect to a point source (§2.2). The photometric behavior near 
a caustic crossing depends on p*, the source size in units of 0 e, whereas the astrometric 
behavior depends on actual angular size of the source 9 * = p*#E. 10 Thus, while fitting to the 
photometric data near a caustic crossing does not yield the angular size of the source, fitting 
to the astrometric data near a caustic crossing does. Therefore it is possible, in principle, 
to measure the angular radius of the source star of a caustic-crossing binary-lens event by 
fitting a few astrometric measurements taken during the caustic crossing to the expressions 
we have derived for the local astrometric behavior (see §2.2.2 and Eq. 60). In practice, 
however, this is complicated by the fact that, for linear fold caustics, only the degenerate 


10 Recall that, in order to covert to astrometric observables, all angular variables must be multiplied by 
0e- 
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combination #*/sin (j) can be measured, where 0 is the angle of the trajectory with respect 
to the caustic. For general fold caustics, the degeneracy between 9 * and sin 0 is broken (see 
§2.2.2 and Eq. 67), however, for many cases the simple linear fold will be applicable. Thus, 
is in order to determine 0* separately, the global geometry must be generally specified, which 
can be accomplished using the global photometric light curve. 


3. A Worked Example: Binary Lensing Event OGLE-1999-BUL-23 

In this section we numerically calculate the expected astrometric behavior for a binary- 
lens fold caustic crossing, and compare this with the analytic results from the previous 
sections. We do this in order to verify our expressions and also to explore the accuracy with 
which our (necessarily) approximate results reproduce the exact behavior. For definiteness, 
we will calculate the expected astrometric behavior for the photometrically well-observed 
caustic-crossing binary lens event OGLE-1999-BUL-23 (Albrow et al. 2001). This has the 
advantage that, up to an orientation on the sky and subject to small errors in the inferred 
parameters, the astrometric behavior can be essentially completely determined from the 
photometric solution, including the size of 9e and the effects of limb-darkening. 


3.1. Formalism and Procedures 

For a system of Ni point masses located at positions Oij, and no external shear, the lens 
equation (Eq. 3) takes the form 


u = 0 


A o-o 

j =i 


id 


1 0-0 


id I 


( 88 ) 


where rrij is the mass of the j'-tli lens in units of the total mass. Note that angles in 
equation (88) are normalized to the 6 e for the total mass of the system. For Ni = 2, the 
lens equation is equivalent to a fifth-order polynomial in 0 , thus yielding a maximum of 
five images. All of the image positions for a given point on the source plane can be found 
numerically using any standard root Ending algorithm. Then the individual magnifications, 
total magnification and centroid of these images can be found using equations (6), (8), and 
(9). 

For a finite source size, it is necessary to integrate over the area of the source (Eq. 49). 
This can be difficult to do numerically in the source place near the caustics, due to the 
divergent magnification. A more robust method is inverse ray shooting. This works as 
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follows. The image plane is sampled uniformly and densely, and at each 6 , equation (88) 
is used to find the corresponding u(0). The local ratio of the density of rays in 6 to the 
density of rays in u is the local magnification. Thus one can create a map of n(u), the 
magnification as a function of u. Similarly, one can determine the astrometric deviation 
by sampling in 0 , using equation (88) to determine u(0), and then summing at each u the 
values of 6 {u). The astrometric deviation at u is then the summed values of 0(u), weighted 
by the local magnification. Thus one creates two astrometric maps, for each direction. In 
practice, inverse ray shooting requires one to bin the rays in the source plane, with the 
resolution of the maps being determined by the size of the bin, and the accuracy determined 
by surface density of rays in the image plane relative to the (unlensed) surface density of rays 
in the source plane. The advantage of inverse ray shooting is that the procedure conserves 
flux, therefore the maps can be convolved with any source profile to produce the finite source 
photometric and astrometric behavior for arbitrary source size and surface brightness profile. 

To apply this method to predict the detailed photometric and astrometric behavior 
of the caustic crossings for OGLE-1999-BUL-23, it is essential that the resolution of the 
astrometric and photometric maps be considerably smaller than the source size p*. We will 
be using resolutions of 1O” 4 0e, which corresponds to 0.034p* for OGLE-1999-BUL-23. This 
is sufficient to accurately resolve the source. We sample the image plane with a density of 
5 x 10% 2 , corresponding to an Poisson error per resolution element of ~ 14%pW Since 
there are ~ 2700 resolution elements per source size, the total Poisson error is always < 1%, 
considerably smaller than any of the effects we will be considering. 


3.2. Global Astrometric Behavior 

Before studying the detailed behavior near the photometrically well-covered (second) 
caustic crossing of OGLE-1999-BUL-23, we first analyze the global astrometric behavior 
of the entire event. We specify the binary-lens topology and source trajectory using the 
parameters of the best-fit solution with limb darkening (see Table 2 of Albrow et al. 2001). 
In Figure 4 we show various aspects of the inferred lensing system and the event itself. The 
best-fit lens system has a wide topology, with two well-separated caustic curves, one near 
the position of each lens. The mass ratio of the two lenses is q = m\jm 2 = 0.39, and they 
are separated by 2.426^. The solution has the source crossing the caustic associated with 
the least massive lens, which we will call the secondary caustic. 

Figure 4(a) shows the photometric light curve centered on the event. There are two 
fold caustic crossings, separated by ~ 10 days; the second crossing was densely covered 
photometrically by Albrow et al. (2001), allowing them to determine not only the source 
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size p*, but also limb-darkening coefficients T in each of two different photometric bands: 
Tj = 0.534 (/-band) and Ty = 0.711 (Id-band). This allows us to predict the astrometric 
behavior including finite source and limb-darkening effects. Furthermore, by combining a 
measurement of p* = with a determination of the angular size of the source from 

its color and magnitude, Albrow et al. (2001) measured the angular Einstein ring radius of 
the lens to be 6e = (634 ± 43)p.as. Therefore we can determine the absolute scale of the 
astrometric features and assess their detectability by comparing them with the expected 
accuracy of upcoming interferometers. 

Figure 4(b) shows the caustics and critical curves, as well as the trajectory of the images 
of the source. One image is always near the most massive lens; this image has little effect 
on the resulting astrometric deviation other than a small net offset along the binary axis. 
Figure 4(c) shows a close-up of the region near the secondary caustic, along with the image 
centroid 6 c i. The components parallel and perpendicular to the binary axis are shown in 
Figures 4(e,d). Finally, Figure 4(f) shows 59 c i, the centroid relative to the unlensed source 
position u. In Figures 4(c-f), the large, 0(6 e) jumps that occur when the source crosses the 
caustic are evident. 


3.3. The Second Fold Crossing 

We now focus on the photometric and astrometric behavior near the second caustic 
crossing for OGLE-1999-BUL-23. In Figure 5(a) we show the caustic geometry and source 
trajectory near this crossing. The source passed ~ 0.02 mas from a cusp. On the scale of 
the source size, 9 * = 1.86/zas, the curvature of the caustic can be neglected. 

In Figure 5 we show the behavior near the second caustic crossing. Figure 5(c) shows 
the photometric behavior for two days centered on the caustic crossing, for the assumption 
of a uniform source, and a limb-darkened source with surface brightness profile given in 
equation (68), with Tj = 0.534 and Ty = 0.711. Note the similarity of the shapes of the 
uniform source and limb-darkened light curves with the analytic forms Go and G 1/2 presented 
in Figure 2(a). We show our prediction for the astrometric behavior as a function of time 
for the components parallel and perpendicular to the binary axis in Figures 5(d) and (e), 
respectively [compare with Figure 1 (e,f)]. We show both the instantaneous discontinuous 
jumps for a point-source centroid, along with the continuous centroid curves for finite sources. 
The limb-darkened and uniform source astrometric curves are extremely similar. In Figure 
5(f) we show the predicted total astrometric behavior for the same time span. Notice how 
well the form of Figure 5(f) compares with Figure 1(g), which was obtained from our analytic 
forms. Figure 6 shows the same curve as Figure 5(f), except we have rotated the axes by 
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~ 55°, shifted the origin to the image of the caustic crossing, and stretched the axes for 
visibility. In Figure 6, notice how the predicted shaped for the finite source centroids are 
smoothly rounded off near the bottom and has a sharp turn near the top. 

Although the difference between the finite source and point source curves is quite sub¬ 
stantial, the difference between the uniform source and limb darkened sources is very small. 
This is more clearly illustrated in Figure 7(a,b), where we show the two components of this 
difference, A 0 l J = 0 l £ — for about four source radius crossing times centered on the 
crossing, i.e. \z\ < 2. The form of AO is very similar to the analytic expectation [see Eq. 86 
and Figure 3(c)]. Furthermore, the two components of AO 1 '} are essentially perfectly (anti-) 
correlated, implying that the difference is essentially one-dimensional. This can be seen best 
in Figure 7(c), where we plot A6 l r '} 2 versus A 0 l ^ 2 . Panel (d) is the same as panel (c), except 
we have rotated and stretched the axes. The maximum deviation is ~ (60r)/xas, in satisfying 
agreement with the rough expectation (Eq. 87). 


4. Summary and Conclusion 

We have presented a detailed study of gravitational lensing near fold catastrophes, 
concentrating on the regime where the individual images are unresolved, i.e. microlensing. 
By Taylor expanding the scalar potential if; in the neighborhood of a fold up to third order in 
the image position, one can obtain a generic form for the lensing map near a fold. Beginning 
with this mapping, we derive the local lensing properties of a source in the vicinity of the 
fold caustic. Approximating the critical curve by its tangent line at the origin, we find that 
the caustic is locally a parabola. On one side of the parabola, the fold lenses a nearby 
source into two images; on the other side of the parabola, there are no images. We derive 
the image positions and magnifications as a function of the position of the source. We find 
that the magnifications of the two images are equal, and recover the well-known result that 
the magnification is inversely proportional to the square root of the distance to the caustic. 
We show how this holds for parabolic caustics (as well as linear caustics), provided that the 
‘vertical’ distance from the caustic is used. 

Assuming a rectilinear source trajectory, and allowing for the existence of slowly- and 
smoothly-varying images not associated with the fold caustic, we derive analytic expressions 
for the total magnification and image centroid (center-of-light) as a function of time. 

We then consider how the photometric and astrometric behavior is altered in the pres¬ 
ence of a finite source size. We derive semi-analytic expressions for the magnification and 
centroid as a function of time for both a uniform source, and limb-darkened source. Along 
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the way we derived expressions that can be used to evaluate the photometric and astrometric 
behavior near a fold for a source with arbitrary surface brightness profile. 

We then show how and under what conditions the generic parabolic fold reduces to the 
more familiar linear fold. We derive simplified expressions for the individual and total image 
positions and magnifications near a linear fold. 

We used some of our analytic results to derive a few generic properties of microlensing 
near folds. In particular, we derive and evaluate expressions for the magnitude of the centroid 
jump that occurs when a finite source crosses a fold relative to the point source jump, and 
the magnitude of the effect of limb darkening on both the photometric and astrometric 
behavior. Notably, we predict, for Galactic bulge lensing events, the shape of the centroid 
due to finite sources with uniform and limb darkening surface brightness profiles. We also 
find for Galactic bulge lensing that the effect of limb darkening on the image centroid near a 
fold is quite similar to the uniform source case, making the limb darkening effect difficult to 
detect by the currently planned accuracy for the instrumentation of SIM. We discussed how 
our formulae can be used to fit both photometric and astrometric data sets near fold caustic 
crossings and thus used to derive such properties as the angular size of the source and the 
microlensing parallax. 

Finally, we numerically calculate expected astrometric behavior of the photometrically 
well-observed Galactic bulge binary lensing event OGLE-1999-BUL-23 (Albrow et al. 2001), 
finding excellent agreement with our analytic predictions. 

Caustics are ubiquitous in gravitational lenses, and the most common type of caustic 
is the fold. Caustics play an especially important role in microlensing, as the rapid time 
variability of the total image magnification allows the possibility of detailed studies of the 
source and lens. In the future, we can expect that time-series photometric measurements will 
be supplemented by time-series astrometric measurements of the center-of-light of microlens 
systems. This paper presents the most thorough and comprehensive study of the photomet¬ 
ric and astrometric behavior of gravitational microlensing near fold caustics to date. The 
results should prove useful to those studying microlens systems with caustics: The analytic 
expressions derived here can be used to fit fold caustic crossings observed both photomet¬ 
rically and astrometrically, gain some insight into more complicated numerical studies, and 
establish predictions for the feasibility of future observations. 
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Fig. 1.— An illustration of the basic properties of astrometric and photometric microlensing 
near folds, (a) Filled circles represent the source at various times. The solid line is the fold. 
The source crosses the caustic at u c . (b) The elongated shapes are the images corresponding 
to the source at the positions in panel (a). The point df c is the image of u Cl and is where 
the two extra images appear. The third image on the nearly horizontal trajectory represents 
the centroid of all images unassociated with the fold, (c) The magnification as a function 
of time for a point source. The solid line is the total magnification g to t, the dotted line is 
the magnification of all images unassociated with the caustic go, and the dashed line is the 
magnification gf of the two images created in the fold crossing, (d) Same as (c), except for 
a finite uniform source. The dashed-dot line shows the total magnification for a source size 
that is two times larger, (e) The 0 c ;,i-component of the centroid shift as a function of time. 
Dotted line is for a point-source, solid line for a finite uniform source size. Dashed-dot line 
is for a source size that is two times larger, (f) The d c ; 2 -component of the centroid shift, (g) 
The path of the centroid of light of all the images, 6 c i. Line types are as in panels (e,f). 
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Fig. 2.— Basic functions which describe the photometric and astrometric behavior of finite 
sources near a fold, as a function of the distance z from the fold in units of the dimensionless 
source size p*. (a) The basic functions for the photometric behavior. The solid line shows 
Go, whereas the dotted line shows G\/ 2 - (b) The basic functions for the astrometric behavior. 
The solid line shows /C 0 , whereas the dotted line shows K, i/ 2 . (c) The photometric (Gi/ 2 — G 0 ) 
and astrometric (/C 1/2 — /C 0 ) limb darkening functions. See text. 
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Fig. 3.— (a) The function [/Co + T(/C i/ 2 — /Co)]/[G*o + r(Gi/ 2 — Go)] as a function of z, 
for T = 0, 0.5, and 1.0, where G 0 , Gn/ 2 , /C 0 , and K-i / 2 are shown in Figure 2, and V is a 
limb-darkening parameter, (b) The fractional difference S/j,^ between the limb-darkened and 
uniform source magnifications, normalized by T, for several source sizes p*. We have assumed 
a caustic scale u/ — 1, and magnification outside the caustic of po,en = 4. (c) The absolute 
magnitude of the difference in the centroid shift due to limb darkening A 0^, normalized by 
TA9j ump , where A 9j ump is magnitude of the point-source astrometric ‘jump’ when the source 
crosses the caustic. See text. 











34 



1300 


1320 


1340 1360 

HJD-2450000. 


1380 


W 

cd 




0 cl l (mas) 


0 cU (mas) 



1300 1350 

HJD—2450000 



Fig. 4.— Global photometric and astrometric behavior of the binary-lens event OGLE-1999- 
BUL-23. (a) The light curve (magnification as a function of time) for the best-fit model. 
Solid line is for the finite source, whereas the dotted line is for a point source. The insets 
show detail near the two caustic crossings, (b) The critical curves (blue ovals), caustics (red 
cuspy curves), images (circles), and source trajectory (dashed line). The X’s denote the 
position of the two masses, the arrows give the directions of motion of the images, and the 
size of the circles are proportional to the magnification of the image, (c) Detail near the 
caustic crossing. The green lines shows the position of the centroid Q c i of the five images 
relative to the lens. (d,e) the two components of 6 c i as a function of time, (f) The solid line 
shows the centroid relative to the source position 50 c i. The X’s show 50 c i at fixed intervals 
of 20 days. 
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Fig. 5.— The photometric and astrometric behavior near the second caustic crossing of 
OGLE-1999-BUL-23. (a) The heavy solid line shows the caustic, while the dashed line 
shows the source trajectory. The interior of the caustic is shaded. The inset shows the detail 
near the caustic crossing, in units of the source size, shown as a circle, (b) The heavy solid 
line is the critical curve, and the circles show the positions of the two images associated with 
the caustic crossing at fixed intervals of 4.8 hours. The size of the circles is proportional to 
the logarithm of the magnification, (c) The magnification near the second caustic crossing 
as a function of time. The solid lines is for a uniform source, dotted line for a limb-darkened 
source in the /-band, dashed line for the E-band, and dashed-dot line is for a point source. 
(d,e) The two components of the centroid Q c i in mas as a function of time. Line types are 
the same as (c). (f) The centroid 6 c i. 
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0’ o u (m as) 

Fig. 6.— Detail of the centroid shift G c i near the second caustic crossing of OGLE-1999- 
BUL-23 in /j, as. This is the same as panel (f) in Figure 5, except the axes have been rotated 
by 55°, and the origin has been translated to the image position of the caustic crossing 
point. Note the extreme asymmetry in the scales of the two axes. The solid line is for a 
uniform source, dotted line for a limb-darkened source in the /-band, dashed line for the 
E-band, and dashed-dot line is for a point source. 
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Fig. 7.— (a,b) The two components of the astrometric offset due to limb darkening A Q 1 ^ 
relative to a uniform source as a function of time. The dotted red line is for the /-band, where 
as the dashed line is for the l/-band. The two components are parallel (a) and perpendicular 
(b) to the binary axis. The vertical lines show (from left to right) the time when the source 
first touches, straddles, and last touches the caustic. The solid lines is the offset from a 
point source, (c) The astrometric offset A O™. Line types are as in panels (a) and (b). 
The deviations from the uniform source are essentially one-dimensional, (d) The same as 
panel (c), except that the axes have been scaled and rotated by ~ 55°. Note the extreme 
asymmetry in the scales of the axes. The maximum absolute deviation due to limb darkening 
is small, < 50/zas. 


















